Comparative study of the rovibrational properties of heteronuclear alkali dimers in 

electric fields 
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A comparative study of the effect of a static homogeneous electric field on the rovibrational 
spectra of several polar dimers in their X^E"*" electronic ground state is performed. Focusing upon 
the rotational ground state within each vibrational band, results for energies and various expectation 
values are presented. For moderate field strengths the electric field-induced energy shifts, orientation, 
alignment, and angular motion hybridization are analyzed up to high vibrational excitations close 
to the dissociation threshold. 

PACS numbers: 33.20.-t,32.60.+i,33.20.Vq 
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INTRODUCTION 

The availability of cold and ultracold dimers provides 
a unique opportunity to investigate fundamental quan- 
tum processes The efforts of many experimental 
groups are focused on the production of heteronuclear al- 
kali dimers, as they constitute a prototype of an ultracold 
quantum gas with long-range dipole-dipole interactions. 
In this context, the photoassociation of LiCs, NaCs, 
KRb, and RbCs has been reported recently 0, H, 0, 0] , 
while Feshbach resonances have been observed for LiNa, 
LiK, LiRb, and KRb (1, 0, H i, [H. The long-range 
anisotropic dipole-dipole interaction between these polar 
molecules gives rise to an intriguing many body physics 
with new interesting physical phenomena. In addition, 
these systems allow for a wide range of applications, such 
as state resolved chemical reactions ll| or quantum com- 
puting 

Particularly interesting is the study of the influence 
of external fields on these molecular systems. Exter- 
nal fields play a key role in trapping and cooling pro- 
cesses as well as in the manipulation and control of the 
dipole-dipole interaction, chemical reactions or collisions 
[14j . As an example, the application of an electric field 
has been predicted to significantly affect the photoasso- 
ciation process of ultracold polar molecules as well as 
the subsequent radiative decay cascade 15, On the 
other hand, the internal structure of molecules is sig- 
nificantly affected by an electric field. The strong field 
regime is characterized by the appearance of pendular 
states: The molecule is oriented along the field axis and 
can be described as a coherent superposition of field- 
free rotational states. Recently, the authors have per- 
formed studies of the rovibrational spectra of CO and 
LiCs in their electronic ground states exposed to an elec- 
tric field: The orientation and hybridization of the angu- 
lar motion, the Stark shift, and the radiative decay prop- 



erties have been investigated for a wide range of field 
strengths and rovibrational levels 17, [3) 
present work we extend our previous studies by perform- 
ing a comparative analysis of several polar alkali dimers 
which are of immediate experimental interest. Specifi- 
cally, we analyze the influence of a static electric field 
on the X^I]+ electronic ground state of ^Li^'^Na, ^Li^"K, 
6Li87Rb, 7Lii33Cs, 23Na39K, 23Nai33Cs, 39K85Rb, and 

85Rbl33Cs. 



ROVIBRATIONAL HAMILTONIAN 

In the following we consider the field regime where per- 
turbation theory holds for the description of the molec- 
ular electronic structure. The Born-Oppenheimer ap- 
proach then provides the following Hamiltonian for the 
nuclear motion of a polar dimer in its electronic ground 
state exposed to a static electric field 



2t2/ 



H = Tn 



:iR) - FD{R) COS 9 (1) 



where the molecule fixed frame with the origin at the 
center of mass of the nuclei has been employed, (i?, 9, (j) ) 
being the internuclear distance and the Euler angles [33|. 
Tfl, ^, hJ{9,(f>), e{R), and D{R) are the vibrational ki- 
netic energy, the reduced mass of the nuclei, the rota- 
tional angular momentum, the field-free electronic po- 
tential energy curve (PEC), and the electronic dipole 
moment function (EDMF), respectively. The electric 
field is taken parallel to the laboratory frame z-axis 
with strength F. Since we restrict ourselves on a non- 
relativistic treatment, couplings to electronic states of 
different spin symmetry do not arise. 

In the framework of the effective rotor approach [l3]. 
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the rovibrational Hamiltonian ([T]) reduces to 
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2n 



[R-^)l3^-F{D{R))lcose + El (2) 



with {R-^)l = mR-^m, {D{R))l = {^^,\D{RM). 

and -E" are the field-free vibrational wave function and 
energy for a given vibrational state v with J = 0, respec- 
tively. This approximation, which is valid for large parts 
of the rovibrational spectrum, relies on two assumptions: 
i) the rotational and vibrational energy scales should dif- 
fer by several orders of magnitude and can therefore be 
separated adiabatically, and ii) the field influence on the 
vibrational motion is very small and can consequently be 
treated by perturbation theory. Besides reducing ([T]) to 
an one-dimensional problem, this approach allows us to 
compare the field influence on several vibrational levels, 
and even between different molecular species, by means 
of the parameter 



h^R 



-2\0 



(3) 



which derives from rescaling the Hamiltonian Q us- 
ing the rotational constant of a given vibrational state, 
B,y — ti^ {R^'^)1/2fi. li u!^ ■ F acquires the same order 
of magnitude as the rescaled rotational kinetic energy, 
, the rotational dynamics is severely altered by the in- 
teraction with the electric field, while aj'^J^ provides an 
estimation of the field strength needed to show a signifi- 
cant effect. 

We remark that in the presence of the electric field, 
only the magnetic quantum number is retained, whereas 
the vibrational and rotational ones are no more con- 
served. However, for reasons of addressability, we will la- 
bel the electrically dressed states by means of their field- 
free vibrational, rotational, and magnetic quantum num- 
bers (vjJjM). Furthermore, in this letter we focus on 
rotationally non-excited states with corresponding field- 
free quantum numbers J = M = 0. The computational 
method used to solve the full rovibrational Schrodinger 
equation is described in ref. [l7j . 



ELECTRONIC POTENTIAL CURVES AND 
DIPOLE MOMENT FUNCTIONS 

The PECs of the singlet electronic ground states of the 
alkali dimers employed here are derived from either ex- 
perimental or theoretical studies. For the PECs of LiNa, 
LiK, LiCs, NaK, NaCs, KRb, and RbCs we use very 
accurate experimental data from spectroscopic analyses 
M, EH, El, lii. El, m, [2^ . The LiRb PEC is taken from 
theoretical studies based on a semiempirical method [13] ■ 
When describing vibrationally highly excited states, spe- 
cial attention has to be payed to the long-range part of 
the PECs: For LiNa, LiCs, NaCs, KRb, and RbCs, the 




FIG. 1: Color online. Electric dipole moment functions of the 
X ^E^ electronic ground states of LiNa (o), LiK (•), LiRb (□), 
LiCs (+), NaK (A), NaCs(v), KRb (*), and RbCs (x). 



PECs include van der Waals terms, — J2n=6 8 10 ^n/R"', 
and an exchange energy term, —AR'^e~^^, see refs. [i^, 

m, Ei, E3, Ea Ei, with 7 



For the LiRb 
-C'e/R^, is con- 



for LiNa [2]|. The 
long-range part of the NaK PEC additionally includes 
damping functions in the van der Waals expansion, 

~J2n=6 8 iQ^n{R)Cn/R^^, and a different expression for 
the exchange energy, —Ae~°'^~^^^ 
dimer only the first van der Waals term, 
sidered |2S |. 

The EDMFs employed for the LiNa, LiK, LiCs, NaK, 
NaCs, and KRb systems are computed by the same 
semiempirical technique as the LiRb PEC [2^. The 
RbCs EDMF is taken from an ab-initio relativistic 
valence-bound calculation [3^. Since for the electronic 
ground state of the polar alkali dimers the long-range 
behaviour of the EDMF is given by Dr/R'^ [3l|, this 
asymptotic function has been fitted to the theoretical 
data at large internuclear distances. We note that, when 
necessary for computational reasons, the short-range be- 
haviour of the EDMFs is linearly extrapolated. The re- 
sulting EDMFs are plotted in fig. [1] as a function of the 
internuclear distance. All EDMFs are negative and ex- 
hibit a broad minimum which is shifted with respect to 
the minimum of the corresponding PEC, e.g., by 0.51 and 
1.39 ao for the KRb and LiCs dimers, respectively. The 
EDMFs of LiRb, NaK, NaCs, and KRb become positive 
for R > 17.3 ao albeit approaching zero thereafter; this 
change of sign will be reflected in the response of highly 
excited rovibrational levels to the field. However, this 
has to be taken with a grain of salt since ab-initio elec- 
tronic structure calculations of molecules typically em- 
ploy basis sets of exponentially localized functions, which 
do not easily catch the long-range behaviour of the atom- 
atom interaction. The LiCs EDMF has the largest ab- 
solute value with a maximum of |Z3(i?)| = 2.36 aoe at 
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FIG. 2: Color online. The parameter |a;^| versus the field- 
free vibrational quantum number v for the same dimers as in 

fig.m 

R = 8.3 ao. The EDMFs of LiNa, KRb, and RbCs show 
only a weak dependence on R for a wide range of inter- 
nuclear distances and overall comparatively small values 
for the dipole moment. 

RESULTS 

Let us start by analyzing the behaviour of the param- 
eter uj^ as a function of the vibrational quantum num- 
ber. In fig. [21 the modulus of uj^, is presented for the 
{ly, 0, 0) levels of all dimers on a logarithmic scale. The 
maximum vibrational quantum numbers considered are 
1/ = 44 (LiK), 46 (LiNa), 48 (LiRb), 53 (LiCs), 73 (NaK), 
85 (NaCs), 97 (KRb), and 135 (RbCs); the very last vi- 
brational bands have not been included. For all dimers, 
\u>^\ shows a qualitatively similar but quantitatively dif- 
ferent behaviour as a function of i^: After a weakly pro- 
nounced maximum, \Lu,y \ decreases rapidly for high vibra- 
tional excitations. In particular, NaK, KRb, and RbCs 
show a plateau-like behaviour indicating a quantitatively 
similar field impact in many vibrational bands. Accord- 
ing to their \uj,y\ values, the considered molecules can be 
classified into four groups whose response to the field will 
be of comparable magnitude for a large part of their spec- 
tra. Possessing the largest \uj^\ values, the first group is 
formed by the NaCs and RbCs dimers and their v < 65 
states, followed by the < 40 levels of LiCs and NaK. 
Despite the large LiCs EDMF, the corresponding \uj^\ is 
only approximately half of the \uj^\ belonging to NaCs 
and RbCs: The large reduced masses of the latter com- 
pensate their smaller dipole moments. The third set con- 
tains LiK, LiRb, and KRb and their lowest 40 vibrational 
bands. Finally, LiNa will be least affected by the field 
which is due to both small D{R) and fi. We remark 



that the sign of Ui, depends crucially on the correspond- 
ing dipole moment: While being generally negative, the 
change of sign of the NaK, NaCs, and KRb EDMFs at 
long distances causes u)^ > for their v > 73, 85 and 95 
levels, respectively. In particular, the minimum observed 
for the V — %5 level of KRb can be attributed to this 
change of sign. For all molecules considered, the slope of 
\uj^ \ becomes very steep as v approaches the dissociation 
threshold. Consequently drastically varies from its 
maximum to its minimum; for example we get for LiNa 
UJ2G = -2.8 X 10"^ versus ujag = -3.0 x 10"^ cm/kV, 
and for NaCs ^^35 = —1.49 versus c^ss = —9.6 x 10"'^ 
cm/kV. Hence, for high-lying vibrational excitations the 
field strength should be enhanced at least one order of 
magnitude to obtain similar field-dressed states as on the 
lower part of the rovibrational spectrum. 

In order to quantify the spectral changes upon the ap- 
plication of the electric field, we consider the ratio of the 
Stark shift to the rotational spacing in the presence of 
the field 

E,,{F)^EMF) ^ ' 

where E^q{F) and E^i{F) are the energies of the {v^ 0, 0) 
and {v, 1, 0) levels at field strength F\ for high-field seek- 
ers we get Ki/ > 0. If <C 1, the field-dressed spectrum 
shows many similarities to the field-free one. In partic- 
ular, this criterion is fulfilled for oj^F ^1: In second 
order perturbation theory the {v, 0, 0) and {v, 1, 0) states 
are shifted by -ujIF'^Bu/Q and [34I, respec- 

tively, compared to the field-free rotational spacing of 

In figs.[3Ja) and (b) the parameter k^, is presented as a 
function of the vibrational quantum number for F = 5.14 
and 51.4 kV/cm, respectively; the selected field strengths 
are chosen such that the experimentally accessible range 
within the strong field regime is covered. Besides show- 
ing high-field seeking character, the behaviour of Ki, has 
many analogies among the considered molecules. Ini- 
tially, it has a weak dependence on the vibrational exci- 
tation, thereby smoothly increasing as v is enhanced and 
passing a shallow maximum (due to the semilogarithmic 
scale, this increasing trend appears as a plateau). For 
higher vibrational excitations, rapidly decreases as the 
levels get closer to the dissociation limit: This reduction 
can be up to one order of magnitude when v is enhanced 
by merely one quantum. The above classification of the 
considered molecules is equally observed for k^, . The low- 
est values appear for LiNa, < 2 x 10"'^ and 0.2, for 
F = 5.14 and 51.4 kV/cm, respectively, i.e., the impact 
of the field on the spectrum is negligible for the smaller 
of these two field strengths. In contrast, NaCs and RbCs 
are most affected by the field F = 5.14 kV/cm, with 
Ki, > 0.5 for 1/ < 68 and 89, respectively. LiCs and NaK 
possess values > 0.2 for < 44 and 58, respectively. 
For LiRb and KRb we found k^, > 0.1 for the < 38 
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FIG. 3: Color online. The relative energy shift 
function of the field-free vibrational quantum number for the 
{v, 0, 0) levels and same molecules as in fig.[T] for (a) F = 5.14 
kV/cm and (b) 51.4 kV/cm. 



and 67 states, whereas a smaller value is achieved for 
LiK with > 0.08 for i' < 34. An increase of the field 
strength by one order of magnitude {F — 51.4 kV/cm) 
produces a significant variation of which depends on 
the system and especially on the states under considera- 
tion. Indeed, we observe > 1 for most levels of LiCs, 
NaK, NaCs, and RbCs. The levels v < 36,41, and 75 of 
LiK, LiRb, and KRb satisfy > 0.8. In order to pro- 
vide a reasonable scaling, the minimum of n^, for KRb at 
v = 95 has not been included in fig. [3l 

The orientation and alignment of the pendular states 
are characterized by the expectation value (cos 9) and 
A cos 6* — (cos^ Q) — (cos^)^, respectively. The closer 
A cos 6* is to zero, the stronger is the alignment and the 
closer the modulus of {cos 9) is to one, the stronger is the 
orientation of the state along the field. The expectation 
value (cos 9) of these levels is presented in figs. IDJa) and 
(b) for F = 5.14 and 51.4 kV/cm, respectively, as a func- 
tion of the field-free vibrational quantum number. The 



corresponding results for Acos0 are shown in figs. [DJa) 
and (b). The sign of the effective electric dipole moment, 
{D{R))1, determines whether the orientation will be par- 
allel or antiparallel to the field direction, i.e., {cos 9) > 
or {cos 9) < 0, respectively. 

For the moment let us focus on the description of the 
results for F — 5.14 kV/cm, see figs. IH^a) and[5l[a). As 
expected from the behaviour of the parameter lo^, the 
orientation smoothly increases as v is enhanced, reaching 
a broad maximum and rapidly decreasing thereafter as 
the state gets close to the dissociation threshold. Since 
the initial increasing trend is not very pronounced for 
NaK, NaCs, KRb, and RbCs, {cos 9) exhibits a plateau 
as a function of v followed by a rapid decrease of the 
orientation. The strongest orientation is achieved for the 
NaCs and RbCs dimers with (cos 61} < -0.7071, for the 
1/ < 65 and 76 levels, respectively. The states v < M 
and 56 of LiCs and NaK, respectively, present a sig- 
nificant orientation as well: {cos 9) < —0.5. For LiK, 
LiRb, and KRb the influence of the field is not very pro- 
nounced and the extremal values are {cos9) = —0.35, 
-0.43, and -0.40 for the v = 22, 23, and 33 levels, re- 
spectively. For LiNa as well as for highly vibrationally ex- 
cited levels of the other molecules, the effect of an electric 
field of strength F = 5.14 kV/cm is negligible. Analo- 
gous conclusions are derived from the analysis of A cos 9, 
whose behaviour as a function of ly is very similar to 
the one previously described for {cos 9). The variation 
of Acos0 is very smooth as v is augmented, and ma- 
jor changes are observed only for high-lying excitations. 
The NaCs and RbCs dimers show the strongest align- 
ment, with A cos 9 < 0.3 for the v <12 and 99 levels, re- 
spectively. The other dimers present relative large values 
with Acos0 > 0.3703 which is the lowest one achieved 
for LiCs at v = 25, indicating a modest alignment. In 
particular, Acos0 is very close to its field- free value, i.e., 
A cos 9 « 0.577, for the LiNa dimer and for all high-lying 
vibrational levels. 

If the field strength is increased to F = 51.4 kV/cm, 
both the orientation and alignment are more pronounced, 
see figs, m^b) and[5^b). Except for the LiNa molecule, 
the strong field regime with the appearance of pendular 
states is already reached for most of the levels consid- 
ered. We observe {cos9) < —0.8 for the states with ly < 
38, 48, 65, 81, and 108 of the LiRb, LiCs, NaK, NaCs, and 
RbCs dimers, respectively. Even more, (cos 6*) < —0.9 
for the 1/ < 70 and 96 states of NaCs and RbCs, re- 
spectively. The LiNa molecule still shows a modest ori- 
entation, with a minimal value {cos9) = —0.394 for the 
ly = 23 state: only if F is further increased by one order of 
magnitude a significant orientation, i.e., {cos9) < —0.8, 
is attained for this molecule. For F = 51.4 kV/cm, the 
strong field regime is also manifest in A cos 9 where many 
levels present a pronounced alignment with A cos 9 < 0.3. 
The smallest values, Acos0 < 0.1, are obtained for the 
1/ < 71 and 100 states of NaCs and RbCs, respectively. In 
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FIG. 4: Color online. The same as fig. |3] but for the expecta- 
tion value (cos 6*) . 



contrast, for the LiNa dimer we find that A cos 6 > 0.485 
for any v. Finally, let us remark that (cos 9) is positive for 
the states with v > 73, 85, and 95 of the NaK, NaCs, and 
KRb dimers, respectively, while being characterized by a 
weak orientation and alignment at both field strengths. 
In particular, for KRb |(cos^?}| (A cos 6') achieves a min- 
imum (maximum) for the ly = 95 level, i.e., for the same 
state as \uj^\ shows a minimum, cf. fig. [2] 

The expectation value (J^) provides an estimate of the 
field-induced hybridization of the angular motion. Since 
we only consider the rotational ground state in this let- 
ter, one yields (J^) = for = 0. Figures El^a) and 
(b) illustrate the evolution of (J^) with increasing degree 
of vibrational excitation for F = 5.14 and 51.4 kV/cm, 
respectively. For both field strengths investigated, (J^) 
presents a similar behaviour for all molecules: It in- 
creases as the vibrational excitation is enhanced, passes 
through a broad maximum, and decreases, approach- 
ing zero, thereafter. The lower parts of the NaK, KRb, 
and RbCs spectra show again a plateau-like behaviour of 
(J^). For F = 5.14 kV/cm, the rotational dynamics has 
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FIG. 5: Color online. The same as fig. |3]but for A cos 6. 



mostly s-wave character, e.g., 47% for the v — 34: level of 
NaCs, 73% for the ly ^ 29 state of NaK, and 87% for the 
J' = 33 level of KRb. The hybridization of the angular 
motion is negligible for LiNa, i.e., (J^) < 3.5 x 10^^ h"^ 
for any v, and 99.8% of the dynamics is controlled by 
the J = wave. In contrast, the pendular character of 
the NaCs and RbCs dimers is already significant with 
(J^) > 1.0 ti^ for the v < 71 and 98 levels, respectively. 
The V < A2 and 52 states of LiCs and NaK, respectively, 
satisfy (J^) > 0.5 /i^; analogously we observe for LiRb 
and KRb (J^) > 0.2 for i/ < 39 and 67. For the LiK 
states we find (J^) > 0.15 for v < 35. 

Because of the high-field seeking character of these lev- 
els, (J^) is drastically increased for = 51.4 kV/cm, see 
fig. \Elh). Except for LiNa, strong field effects charac- 
terize this expectation value for a wide range of levels. 
Indeed, (J^) > 5.0 /i^ for the v < M and 73 states of 
NaCs and RbCs, (J^) > 3.0 fi^ for those with v < 50 
and 41 of the NaK and LiCs dimers, and (J^) > 1.8 ft^ 
for the lowest 35 vibrational states of LiK. They present 
a rich pendular dynamics, with significant contributions 
of higher field-free rotational states. As v gets closer to 
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FIG. 6: Color online. The same as fig. |3] but for the expecta- 
tion value (J^). 



the dissociation limit the relative weight of the s-wave 
becomes more important, and (J^) is reduced. In con- 
trast, the field- free J = level still dominates the LiNa 
states with (J^) < 0.26 for any z/- value. Finally, let us 
remark that the effect of the field in terms of (J^) for the 
v = 95 level of KRb is negligible but slightly increases 
again for the last states. 



CONCLUSION 

In conclusion, we have investigated the rovibrational 
spectra of the electronic ground state X^E+ of the al- 
kah dimers LiNa, LiK, LiRb, LiCs, NaK, NaCs, KRb, 
and RbCs in the presence of a static homogeneous elec- 
tric field. Performing a comparative analysis, we have 
considered the lowest rotational excitation within each 
vibrational band up to the dissociation limit. The vi- 
brational state-dependent parameter Ui, provides a first 
estimation of the influence of the field on the correspond- 
ing dimer and allows a classification of these systems into 



four groups. Thus, NaCs and RbCs are most affected by 
the application of an electric field and LiNa is least. In 
general, vibrationally highly excited states show less field 
impact whereas already for moderate field strengths pro- 
nounced effects are observed for lower-lying levels. The 
richness and variety of field-dressed rotational dynam- 
ics have been illustrated by analyzing the Stark shift of 
the spectra as well as the orientation, alignment, and 
hybridization of the angular motion. 

This work should serve as a guide to experimentalists 
interested in electric field effects on cold and ultracold 
heteronuclear alkali dimers. A natural extension of this 
work includes the study of the influence of the field on 
states emerging from higher rotational excitations and/or 
with nonzero magnetic quantum numbers M. 
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